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Highly efficient broadband polarization retarders and tunable polarization filters 
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By using the formal analogy between the evolution of the state vector in quantum mechanics 
and the Jones vector in polarization optics, we construct and demonstrate experimentally efficient 
broadband half-wave polarization retarders and tunable narrowband polarization filters. Both the 
broadband retarders and the filters are constructed by the same set of stacked standard multi-order 
optical wave plates rotated at different angles with respect to their fast polarization axes: for a 
certain set of angles this device behaves as a broadband polarization retarder while for another 
set of angles it turns into a narrowband polarization filter. We demonstrate that the transmission 
profile of our filter can be centered around any desired wavelength in a certain vicinity of the design 
wavelength of the wave plates solely by selecting appropriate rotation angles. 

PACS numbers: 


I. INTRODUCTION 

Optical polarization retarders have countless appli¬ 
cations in physical experiments. Retarders are opti¬ 
cal plates, made of birefringent material, which possess 
two orthogonal axes — ordinary and extraordinary — 
with different refractive indices, rig and rig. Light of 
wavelength A passed orthogonally through a retarder of 
thickness L experiences a retardation (phase shift) ip = 
2TTL{ne—no) /A. Common types of wave plates (WPs) are 
half-wave plates (HWP) with retardation (p = {m + l)7r 
and quarter-wave plates (QWP) with (p = (m -I- ^)7r, 
(m = 0,1, 2,...). The QWPs and HWPs can be zero- 
order (m = 0) or multi-order (m > 0). Because p de¬ 
pends on A, the retarders exhibit wavelength dispersion 
and therefore QWP and HWP are produced with a spec¬ 
ified design wavelength. Application^ e.g. in the ter¬ 
ahertz time-domain spectroscopy [I|, 3, microwave po- 
larimetry etc. often demand wider wavelength 

range of nearly constant phase retardation. 

Achromatic (or broadband, BB) retarders were first 
proposed by West and Makas Q by combining two or 
more plates with different birefringence, whereby one of 
the plates is rotated to a specified angle. Destriau and 
Prouteau [3] used two WPs of the same material but 
with different thickness and axis angles. Pancharatnam 
used three plates to construct half-wave @ and quarter- 
wave 0 BB retarders. Harris et al. proposed achromatic 
QWPs with 6 and 10 [ll| stacked identical zero-order 
QWPs rotated at different angles, optimized to reduce 
the wavelength dispersion over an extended range. Alter¬ 
natively, for some applications tunable phase retardation 
plates can be used. Polarization retarders can be made 
tunable either by tilting the WP by a specified angle (Al- 
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phalas GmbH, Gottingen) or electrically, by changing the 
birefringent optical path length using a twisted nematic 
liquid crystal phase retarder Bill- 

Birefringent phase retarders have been also used to 
construct narrowband (NB) optical filters. The Lyot 
BB and Sole filters B, B, B build up by stacks 
of retarders and polarizers. Possibilities to obtain very 
narrowband transmission filters for applications, e.g. so¬ 
lar ima ging at specific wavelengths, were also demon¬ 
strated [isi B ■ 

Electrically tunable polarization optical filters based 
on the Lyot filter scheme, using the birefringent effect 
of electrically controlled liquid crystalline la yers instead 
of solid retarders have been investigated [20l - B and are 
already commercially available. Due to the use of polar¬ 
izers and the limited transmission of liquid crystals, how¬ 
ever, the overall peak transmission of filters using higher 
number of polarizers and/or liquid crystal elements is 
relatively poor. 

Stacked composite plates are mathematically equiva¬ 
lent to composite pulses in quantum physics, which are 
widely used in nuclear magnetic resonance (NMR) B 
and since recently, in quantum optics IB - This 
similarity stems from the formal analogy between the 
Schrodinger equation for a two-state quantum system 
and the equation for the evolution of the Jones polar¬ 
ization vector B- If Ii9-s been used by Ardavan B 
who proposed to use the BB and NB composite pulses of 
Wimperis B to design BB and NB WPs. In a recent 
work (B - ths analogy between the polarization Jones 
vector and the quantum state vector was used to pro¬ 
pose a calculation method to derive the twisting angles 
of a stack of retarders with arbitrarily accurate BB po¬ 
larization retarders, which promise to deliver very high 
polarization conversion fidelity in an arbitrarily broad 
range of wavelengths. 

In this paper, we follow the theoretical proposal of 
Ivanov et al. [32j| to design and demonstrate experimen- 
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tally stacked composite plates, which, when rotated at 
specific angles, act either as BB half-wave retarders or 
tunable NB filters. In contrast to the earlier demonstra¬ 
tion by Peters et al. [s^, here we use the more flexible 
single-pass setup rather than the double-pass setup used 
there. Moreover, we use multiple-order WPs, which are 
cheaper and easier to produce, and are therefore advan¬ 
tageous for practical reasons. We also use a lamp as a 
light source, rather than lasers as used in Ref. [s^ , which 
allows us to obtain more detailed transmission profiles 
in a wider wavelength range. Furthermore, the devel¬ 
oped model shows how to tune the NB filter to any de¬ 
sired wavelength (other than the design wavelength Aq) 
by merely selecting suitable rotation angles. 

In Sec. m we present the theoretical method used to 
design our polarization BB retarders and NB filters. The 
experimental apparatus is described in details in Sec. imi 
In Sec. lIVi we show the measured transmittance spectra 
of our composite filters and BB retarders. In this section, 
we also show how to tune the central wavelength of the 
NB filters by simply selecting suitable rotation angles. In 
Sec. El we present the conclusions. 


II. THEORY AND NUMERICAL 
CALCULATIONS 


Here we will summarize the basic theory of compos¬ 
ite optical retarders. In the Jones calculus [2^, a single 
birefringent retarder is described by the matrix 


d a (^) 


R{—a) 


0 

0 


R{a), 


where 


R{a) 


cos a sin a 
— sin a cos a 


(1) 

( 2 ) 


We use the linear polarization basis, a pair of orthogonal 
polarization vectors, where a is the rotation angle of the 
retarder’s optical axis and ip = (/^(A) is the retardation 
accumulated between the ordinary and the extraordinary 
rays passing through the retarder. A stack of iV retarders 
is described by the composite Jones matrix 


((^^) (</57V-i) ■ • ■ Jai (V^l) , (3) 


where the light passes through the plates described by 
Jqj Wi) in ths order of ascending index j, i.e. from right 
to left, and $ is the overall composite retardation. For 
the sake of brevity, we denote ($) = J. For practical 
purposes we take all plates to be QWP (Q) with respect 
to the design wavelength Aq . Thus we set pj = po , where 
‘Po = 7!‘/2 for Aq. This allows us to use commercially 
available standard WPs. 

We thereby assume the following configuration of N 
plates each rotated at an angle aj (cf. Fig. [T|), 


(4) 


The overall composite retardation is $ = 2 cos“^ Re Jn. 
The element Jn shows the fraction of light intensity 
which survives when light of particular polarization is 
passed through the filter. The retardation profile is sym¬ 
metric relative to the target retardation $ = tt, i.e., 
$((/?) = $(27r — p). 
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FIG. 1; Experimental setup. The source S, irises 7i and I 2 , 
lens Li and polarizer Pi form a collimated beam of white 
polarized light. Polarizer P 2 and lens L 2 focus the beam of 
output light onto the entrance F of an optical fibre connected 
to a spectrometer. The composite retarders and filters are 
constructed by a stack of multiple-order QWPs (Qn). 

By varying the JV rotation angles aj we obtain different 
half-wave composite BB retarders and NB filters. 

For RB retarders the angles aj satisfy 

max(|$((/j)/7r - 1|^) < e, p & [puiir^T^ - Pmm] ■ (5) 

This guarantees that there is a range of single plate 
retardations between i^min and tt — pmin (with p = tt 12 
corresponding to Aq), where the composite retardation 
$ remains close to tt. We determine numerically those 
angles aj, which minimize Pmin, giving the broadest pos¬ 
sible range of high-quality retardation. 

On the contrary, a NB filter must transmit only a small 
spectral region around Aq. This imposes the following 
relation on the angles Uj: 

max(| Ji 2 |^) < e, lyS G [0, U [tt - (/Jmax, tt] , (6) 

where |Ji 2 |^ represents the filter’s intensity transmit¬ 
tance from horizontal to vertical polarization, as known 
from Jones calculus. We seek those angles aj, which 
maximize v^max, thereby giving the narrowest possible 
transmission window around 7r/2. Outside this window 
the transmission is maintained below e. For both the BB 
and NB retarders we choose e = 1%. 

We set the optical axes of the composite retarders and 
filters at angle 0° by setting |Ji 2 | = 1 a,i p = 7r/2. 
This imposes an additional constraint on the angles aj 
(cf. Fig. [T]), thereby leaving N — 1 angles aj free to vary. 

The more angles we can vary, the easier it is to ful¬ 
fill Eqs. ([S]) and ([5]). Therefore longer sequences (|1]) 
of larger number N of constituent wave plates provide 
larger bandwidths of the broadband retarders or smaller 
bandwidths for the narrowband filters. 

It is remarkable that we can tune the wavelength of 
the filter, i.e. we can center the transmission window 
at a wavelength A' other than the design wavelength of 
the single plate Aq, just by varying the values of the ro¬ 
tation angles a^. The latter are obtained from Eq. ©, 


Qa,Qa.- -Q 
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where now an offset in the argument of J 12 is introduced: 

t Here ^p' corresponds to the desired 

new wavelength A': Lp' = 2TrL{ne — no)/X'■ 

For the numerical optimization we use Newton’s 
gradient-based method. Numerous solutions to Eqs. 
or (|6]) exist depending on the values of v^min or ipmax- To 
determine the optimal solutions we gradually decrease 
P>min or increase until we reach the least or the 

largest ipmax, which still satisfy Eqs. m or ®. Because 
we use a local optimization algorithm, we iteratively pick 
the initial parameter values using a Monte-Carlo scheme. 
Although the calculations are made for zero-order plates 
(m=0), they are valid for multi-order plates too as the 
matrix Jai^p) and the ensuing matrix are peri¬ 

odical in the retardation ip. 

The calculated results for the rotation angles aj of the 
optical axes of the individual plates for several compos¬ 
ite sequences are given in the captions of the respective 
figures below. 


III. EXPERIMENT 

We have tested experimentally the calculated BB re¬ 
tarders and NB filters by using the setup shown in Eig.[T] 
A 10 W Halogen-Bellaphot (Osram) lamp with DC power 
supply was used as a source of white light. A collimated 
light beam was made by using a variable iris Ii with 
an aperture less than 0.5 mm positioned at the focus of 
a plano-convex lens Li (/=150 mm). A second iris di¬ 
aphragm I 2 was used to avoid outer rings of the beam 
pattern. The diameter of the beam, measured at a dis¬ 
tance of 2 m, was about 2 mm. The light beam was po¬ 
larized in the horizontal plane by a linear polarizer Pi. A 
second polarizer P 2 was placed at a distance of 0.7 m and 
served as an analyzer during the experiments. Both po¬ 
larizers (Clan-Tayler, 210-1100 nm) were borrowed from 
a Lambda-950 spectrometer. A second plano-convex lens 
L 2 (/=20 mm) and a two-axis micro-positioner were used 
to focus the light beam onto the optical fibre entrance F 
connected to a grating monochromator (Model AvaSpec- 
2048 with controlling software AvaSoft 7.5). Multi-order 
QWPs (WPMQlOM-780, Thorlabs) were used to build 
the composite retarders as described above. At wave¬ 
length 780 nm the multi-order QWPs are designed to 
be 11.25 waves, while at 765 nm they work as HWPs. 
Each WP (aperture of 1") was assembled onto a sepa¬ 
rate RSPl rotation mount. The optical axes of all the 
WPs were determined with an accuracy of 1°. With the 
used light source and monochromator, we could obtain 
reliable spectral transmission data in the range of 400-830 
nm. 

The composite retarders were assembled as a set of 
multi-order WPs, each of them rotated at specific angles 
provided by the calculations, as described above. The 
WPs were slightly tilted [1^ to reduce unwanted reflec¬ 
tions. We used a single beam spectrometer, thus all ex¬ 
periments started with measurement of the dark and ref¬ 
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FIG. 2: Measured transmittance for BB composite HWP 
formed of a sequence of QWPs: a) reference spectrum of a 
HWP made of two QWPs; b) BB HWP retarder made of 6 
QWPs with rotation angles (45.5°, 78.5°, 76.7°, 15.5°, 17.7°, 
45.4°). 


erence spectra. The dark spectrum, taken with the light 
path blocked is further automatically used to correct for 
hardware offsets. The reference spectrum is usually taken 
with the light source on and a blank sample instead of the 
sample under test. In our case, however, we have mea¬ 
sured the transmission spectrum of the already assembled 
composite retarder, but the axes of the polarizer Pi, the 
analyzer P 2 and the fast axis of the single WPs were all 
set parallel. The measured light spectrum was used as 
a reference for the subsequent measurements with the 
waveplates rotated at their respective theoretically cal¬ 
culated angles aj, and the analyzer P 2 set to 90°. The 
transmission spectrum of the so-obtained composite re¬ 
tarder was recorded and scaled to the reference spectrum. 
The unavoidable losses due to reflections and absorptions 
from any single waveplate were thereby taken into ac¬ 
count. The obtained spectrum predominantly depends 
on the effective retardance of the composite retarder. 

When the composite half-wave retarder was built by 
QWPs, the measured transmission spectra for both BB 
and NB retarders were compared with the spectrum of 
two QWPs assembled to act as a single HWP. In the cases 
when the WPs WPMQlOM-780 were used as HWPs at 
765 nm, the measured transmission spectrum of a single 
WP was used to demonstrate the broadening or narrow¬ 
ing effect. 

According to the theoretical model there are many 
suitable sets of angles for the composite retarders. In 
all cases we have tested experimentally up to 10 differ¬ 
ent sets of angles for each retarder type, but here, for the 
sake of brevity we present only the representative results. 
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FIG. 3: Measured transmittance for BB composite HWP as¬ 
sembled by a sequence of QWPs used as HWPs at 765 nm. 
a) reference spectrum of one QWP; b) BB HWP retarder 
made of 3 QWPs with rotation angles (14.7°, 164.4°, 14.7°); 
c) BB HWP retarder made of 5 QWPs with rotation angles 
(7.5°, 172.5°, 14.2°, 172.9°, 8.6°). The inset presents the re- 
tardance. 


IV. EXPERIMENTAL RESULTS 
A. Broadband polarization retarders 

In this section we will demonstrate that using the cal¬ 
culation method described above composite half-wave 
retarders with enhanced wavelength range can be as¬ 
sembled using a set of multi-order WPs. The design 
wavelength of the used multi-order (11.25 waves) QWPs 
(WPMQlOM-780) is 780 nm. These WPs could also 
be used at 765 nm as HWPs (cf. the data sheet of 
WPMQlOM-780, Thorlabs). We will demonstrate here 
how these commercial QWPs can be used to construct 
a composite BB HWP for both the designed wavelength 
780 nm and 765 nm. 

First, a composite HWP was made using six multi¬ 
order QWPs WPMQlOM-780. The QWPs were twisted 
to the respective angles listed in the caption of Fig. [21 
Figure 12] shows the transmittance curve. For comparison, 
the transmittance spectrum of a HWP build up of two 
QWPs with parallel axes is shown too. The significant 
wavelength range broadening yielded with the composite 
HWP retarder is evident: the width of the 90%-efficiency 
range is nearly tripled. 

In a second experiment, a BB HWP for 765 nm was 
composed using the same QWPs as above, but now used 
as HWPs at 765 nm. Representative transmittance spec¬ 
tra of two composite BB HWPs, assembled with 3 and 5 
HWPs, respectively, are presented in Fig. [31 The trans¬ 
mittance spectrum of a single HWP at 765 nm is shown 
for comparison. The calculated retardance is shown as 
inset in Fig. 121 As in Fig. |21 we find significant wave¬ 
length range broadening. As expected, the composite 
HWP with higher number of constituent WPs has a wider 
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FIG. 4: Measured transmittance for NB composite filters, 
a) reference spectrnm of two QWPs; b) filter of 5 QWPs, 
rotation angles (43.7°, 176.9°, 170.1°, 119.3°, 80.2°); c) filter 
of 6 QWPs, rotation angles (165.3°, 167°, 19.7°, 18.6°, 166.4°, 
166.1°). 


wavelength range. 


B. Composite polarization filters 

Standard multi-order WPs are NB retarders. Using a 
proper set of rotation angles, a composite NB retarder 
can be assembled with much narrower wavelength range. 
Here we demonstrate how with the same set of QWPs, 
as the one used above, but rotated at different angles, 
a NB HWP can be assembled. Figure [4] shows repre¬ 
sentative transmittance curves for composite NB filters 
constructed of 5 and 6 QWPs, respectively. The 6-plate 
composite NB filter reduces the width (at half-maximum) 
of the transmittance spectrum by a factor of almost 3. 
The comparison of the transmittance spectra (curves b 
and c) shows that the filter bandwidth is narrower for a 
higher number of WPs. Even with this limited number of 
WPs one can produce filters with a bandwidth of about 
7 nm. There is no fundamental physical limit (except the 
divergence of the light beam), which can prevent us from 
going to as narrow width of the transmission profile as 
we like. With sufficiently many WPs, one can reduce the 
bandwidth even below 1 nm for sufficiently collimated 
light beam. 

As for BB retarders, we have calculated sets of angles 
for the NB HWPs assembled using the QWPs as HWPs 
at 765 nm. The results are shown in Fig. |5l Again, a 
reduction of the transmittance bandwidth by a factor of 
3 is established with the 6-plate composite filter. 
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FIG. 5: Measured transmittance for composite filters, a) ref¬ 
erence spectrum of one HWPs; b) filter of 3 HWPs, rotation 
angles (14.7°, 45.1°, 75.3°); c) filter of 5 HWPs, rotation an¬ 
gles (37.7°, 83.1°, 49.9°, 1.9°, 38.2°). 


C. Tunable composite polarization filters 

As already shown, very narrow-band optical filters can 
be assembled using a stack of QWPs or HWPs with prop¬ 
erly chosen rotation angles. Furthermore, the calculation 
method presented above can be also applied to solve the 
problem of making these filters tunable. The tunabil- 
ity was tested for a composite HWP filter composed of 
6 QWPs. The measured set of transmittance spectra is 
shown in Fig. (5] It is seen that by choosing properly the 
rotation angles the filter’s spectral band can be moved 
at will. Remarkably, the bandwidth of the transmittance 
curves does not change with the tuning of the filter. Even 
more, one could use the multi-order WPs as QWPs or as 
HWPs at different wavelengths as shown in the previous 
section, thus constructing tunable narrow-band filters in 
a wide wavelength range. 


at 765 nm and QWP at 780 nm, one can build effi¬ 
cient polarization BB (achromatic) half-wave composite 
retarders and tunable NB polarization filters. By rotat¬ 
ing the same set of WPs at specific angles we have con¬ 
structed BB HWP around 780 nm, BB HWP around 765 
nm, and tunable HWP filter at several wavelengths from 
765 nm to 788 nm. Better (broader) BB retarders can 
be built with lower-order WPs. Polarization filters with 
narrower bandwidth (FWHM), even below 1 nm, can be 
constructed with sufficiently many high-order (the higher 

the better) WPs. 
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FIG. 6: Measured transmittance of composite filters made of 
the same set of 6 QWPs but for different rotation angles. The 
central wavelength is tuned to: a) 780 nm (design wavelength, 
no tuning, retardation ip' = 0.507r); b) 772 nm [p' = 0.257r); 
c) 775 nm {p' = 0.357r); d) 784 nm {p' — 0.657r); e) 788 nm 
{p' = 0.757r). The rotation angles are: a) (165.3°, 167.1°, 
19.7°, 18.6°, 166.4°, 166.1°); b) (8.9°, 158.1°, 158.9°, 116.1°, 
51.5°, 34.5°); c) (26.4°, 154.7°, 178.6°, 1.3°, 20.2°, 158.8°); d) 
(60.4°, 14.1°, 175.7°, 178.4°, 154.2°, 110.6°); e) (6.2°, 25.4°, 
50.2°, 56.1°, 39.0°, 13.4°). 
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